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[APPLICANT HEREBY REQUESTS PARTICIPATION IN THE PATENT PROSECUTION HIGHWAY {PPH} pilot program 
AND PETITIONS TO MAKE THE ABOVE-IDENTIFIED APPLICATION SPECIAL UNDER THE PPH PILOT PROGRAM. 



The above-identified application (1) validly claims priority under 35 U.S.C. 119(a) and 37 CFR 1.55 to one or more 
corresponding IL appllcation(s) or to a PCT application that does not contain any priority claim, or (2) is a national stage 
entry of a PCT application that does not contain any priority claim. 

The IUPCT application n */nnnnft 
numbers) is/are: PCT1L05/Q001 lb 

IL 160152 



The filing date of the IU ^ _„ _ 4 or . nf ~ 

pct applications) istaw PCT filing date- February 1 ,2005 

IL filing date- February 1 , 2004 

i. List of Required Documents; 

a, A copy of alt IL office actions which are relevant to patentability in the aboveidentifled IL application^} 



is attached. 



b A copy of all claims which were determined to be patentable hy the IL in the above-identified IL 
application^) 

\& attached. 

c English translations of the documents in a. and b. above along with a statement that the English 
translations are accurate are attached (if the documents are not in the English language). 



d. (1) An information disclosure statement listing the documents cited in the IL office actions 

□ Is attached- 

ly Has already been Hied in the above-identified U.S. application on 10/31/2010 



{2) Copies of ail documents {except for U.S. patents or U.S. patent application publications) 

□ Are attached. 

✓ t Have already been filed in the above-identified U.S. application on 1 0/31/ 201 0 
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Hi. All the claims in the US application sufficiently correspond to the patentabie/allowable clatrns in the IL 
application. 
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Privacy Act Statement 



The Privacy Act of 1974 (PX. 93-579) requires that you be given certain information in connection with your 
submission of the attached form related to a patent application or patent. Accordingly, pursuant to foe requirements of 
the Act T please be advised that; (1) the general authority for the collection of this information rs 35 U.S.C. 2(b)(2); (2) 
furnishing of the Information solicited is voluntary; and (3) the principal purpose for which the information is used by the 
US. Patent and Trademark Office is to process and/or examine your submission related to a patent application or 
patent, if you do not furnish the requested information, the U.S. Patent and Trademark Office may not be able to 
process and/or examine your submission, which may result in termination of proceedings or abandonment of the 
application or expiration of the patent. 

The information provided by you in this form wiB be subject to the following routine uses: 

1 . The information on this form wiii be treated confidentially to the extent allowed under the Freedom of 
Information Act (5 U.S.C, 552) and the Privacy Act (5 U.S.C 552a). Records from this system of records 
may be disclosed to the Department of Justice to determine whether disclosure of these records is required 
by the Freedom of Information Act, 

2. A record from this system of records may be disclosed t as a routine use, tn the course of presenting 
evidence to a court magistrate, or administrative tribunal, including disclosures to opposing counsel in the 
course of settlement negotiations, 

3. A record in this system of records may be disclosed, as a routine use, to a Member of Congress submitting 
a request involving an individual to whom the record pertains, when the individual has requested 
assistance from the Member with respect to the subject matter of the record. 

4. A record in this system of records may be disclosed, as a routine use, to a contractor of the Agency having 
need for the information in order to perform a contract Recipients of information shall be required to 
comply with the requirements of the Privacy Act of 1974, as amended, pursuant to 5 U.S.C, 552a(m). 

5. A record related to an International Application fried under.the Patent Cooperation Treaty in this system of 
records may be disclosed, as a routine use, to the Internationa! Bureau of the World Intellectual Property 
Organization, pursuant to the Patent Cooperation Treaty. 

6. A record in this system of records may be disclosed, as a routine use, to another federal agency for 
purposes of National Security review (35 U.S.C. 181) and for review pursuant to the Atomic Energy Act (42 
U.S.C. 218(c)). 

7. A record from this system of records may be disclosed, as a routine use, to the Administrator, General 
Services, or his/her designee, during an inspection of records conducted by GSA as part of that agency's 
responsibility to recommend improvements in records management practices and programs, under 
authority of 44 U.S.C. 2904 and 2906. Such disclosure shaH be made in accordance with the GSA 
regulations governing inspection of records for this purpose, and any other relevant ( i.e. , GSA or 
Commerce) directive. Such disclosure shall not be used to make determinations about individuals. 

8. A record from this system of records may be disclosed, as a routine use, to the public after either 
publication of the application pursuant to 35 U.S.C, 122(b) or issuance of a patent pursuant to 35 U.S.C. 

1 51 , Further, a record may be disclosed, subject to the limitations of 37 CFR 1 .14, as a routine use* to the 
public if the record was filed in an application which became abandoned or in which the proceedings were 
terminated and which application is referenced by either a published application, an application open to 
public inspection or an issued patent. 

9. A record from this system of records may be disclosed, as a routine use, to a Federal, State, or local law 
enforcement agency, if the USPTO becomes aware of a violation or potential violation of law or regulation. 
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Israeli Patent Office 



22 June 2008 



Notification of deficiencies in Israeli application 160152 

In correspondence with rule 41 of the Patent Rules, we hereby inform you that deficiencies 
have been found in the above-mentioned application. 

A response must be submitted within 4 months, extensions of time requested before the 4 
months' deadline at the cost of 56 NIS per month. 

1. The set of claims includes too many independent claims for method (4), in 
contradiction to what is customaty, see Registrar decision in opposition 1 12858 
dated November 29 2004, clauses 35-38 ( and also see rule 20 (a) (3)).Allegedly they 
refer to distinct inventions in contradiction to clause 8 of the Law, 

2. In order to save time, a preliminary search has been carried out, Claim 1 and 8 lack 
inventive step, as required by clause 5 of the Patent Law, in light of the following 
publication: REICHENBACH "Restoration and Reconstructions from Overlapping Images 
for Multi-Image Fusion" IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE 
SENSING< VOL 39, NO 4 April 2001, pp 769-780 , which describes a method for 
reconstructing a scene from multiple samples, which includes a weighted average of 
values corresponding to the spatial location in relation to the output ( see page 778, right 
column, line 26 until page 779 left column, line 2 and figure 12). Combining this approach 
with the process of super resolution (which is a particular example of such reconstruction) 
is described in the prior art publications mentioned in the summary, and are obvious to a 
person skilled in the art. 

The same applies, with the relevant changes, to claim 18, which refers to a similar 
combination, 

3. Please submit an update of prior art references cited in corresponding 
applications. 

Respectfully, 
Eran Ross 
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IN THE ISRAEL PATENT OFFICE 

26 May 2009 
Our File: 26/604 

Memorandum re: Super-Resolution Image 
Processing 

Israel Application No A 60 1 52 

Filed: 1 February 2004 

By; Rafael Advanced Defense Systems 

Ltd* 

Dear Sirs, 

This is a response to the Examiner's letter of 22 June 2008, which is being 
filed on or before 28 May 2009 and for which eight months' extension fees have been 
paid. 

Enclosed please fmd replacement pages 22-27 of the specification, in 
duplicate, containing claims 1-11 which replace pages 22-32, claims 1-25, as filed. 

The Examiner has objected to the presentation of 4 independent method 
claims, citing the Patent Registrar's decision in the opposition to Israel Patent 
Application 1 12858. 

In order to address this issue, the Applicant has chosen to cancel claims 1-6 
and 18-25 from the application, without prejudice. The remaining claims 7-17 are 
unchanged, but have been renumbered as claims 1-11. The Applicant reserves the 
right to file a divisional application under Section 24(h) of the Israel Patent Law 
(1967), 

The Examiner has provisionally raised a question under Section 8 of the Law 
as to whether the claims relate to more than one invention. In response, the Applicant 
respectfully submits that the independent claims originally presented relate to a 
number of closely interrelated aspects which are preferably used together in synergy. 



i 
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Thus, for example, the distinguishing features of original independent claim 13 (now 
claim 7) also appear in dependent claim 10 (now claim 4) depending from original 
independent claim 7 (now claim 1), while the distinguishing features of original claim 
7 also appear in original dependent claim 16 (now claim 10) depending from original 
claim 13 (now 7). 

Given that the various features recited in the independent claims all relate to 
the same overall process, and recite closely interrelated features which are preferably- 
used together in synergy, the Applicant respectfully submits that the claims should not 
be considered "different inventions" in the sense of Section 8 of the Israel Patent Law. 

The Examiner has also cited a Reichenbach reference as relevant to original 
claims 1 and 18- The Applicant respectfully disagrees. However, this issue is no 
longer relevant in view of the cancellation of claims 1-6 and 18-25 from the 
application. 

Finally, regarding the Examiner's request under Section 18, we have received 
an International Search Report in the parallel PCT application no. PCT/ILQ5/000116. 
A copy of the search report and an accompanying CD are being submitted as part of 
an updated Section 18 submission in parallel to this response. 

In view of the accompanying amendments and remarks, the Applicant believes 
that all issues raised by the Examiner have been folly addressed. Positive 
consideration is respectfully solicited. 

Respectfully submitted, 

Dr. Mark M. Friedman 
Advocate, Patent Attorney 
Mo she Aviv Towers, 54 lh Floor 
7 Jabotinski Street 
Ramat Gan 52520 
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WHAT IS CLAIMED IS: 

1. A method for iterative derivation of a master image from a 
plurality of sampled images of non~ identical, at least partially overlapping, 
regions of a scene, the master image having an output resolution greater than a 
maximum resolution of the sampled images, the method comprising; 

(a) for each sampled image, defining: 

(i) a transformation operator F mapping positions within the 
master image to corresponding positions in the sampled 
image, 

(ii) a distortion operator H simulating a distortion associated 
with an imaging sensor ftom which the sampled image 
was generated,, and 

(iii) a sampling operator D for reducing an image from the 
output resolution to the resolution of the sampled image; 

(b) for each sampled image, applying said transformation operator, 
said distortion operator and said sampling operator to a current 
master image hypothesis so as to generate a predicted image, and 
calculating a difference image having pixel values corresponding 
to the difference in corresponding pixel values between the 
sampled image arid the predicted image; 

(c) performing back-projection of each of said difference images to 
generate a correction image for the current master image 
hypothesis; and 
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(d) employing said correction images to perform a coxrection to the 
current master image hypothesis to generate a new master image 
hypothesis, 

wherein said back projection includes employing m operator H bp 
corresponding to a pseudo-inverse of distortion operator H, wherein H bp 
approximates to an inverse of H at spatial frequencies below a given value and 
approaches zero at spatial frequencies above said given value. 

2. The method of claim 1, wherein H bp is chosen to substantially 
satisfy the condition; 

wherein: 

I is the unit operator for an image of the output resolution; 

D is a sampling operator for reducing an image from the output 

resolution to the resolution of an input image; and 
D is an inflation operator for expanding an image from the resolution of 

the input image to the output resolution. 

3. The method of claim 1, wherein distortion operator H 
corresponds to a combination of a modulation transfer function resulting from 
an optical system of the imaging sensor and a modulation transfer function 
resulting from a distortion generated by a sensor element array of the imaging 
sensor. 
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4. The method of claim 1, wherein distortion operator H 
corresponds to a modulation transfer function describing only a first portion of 
a distortion generated by the imaging sensor, the method further comprising a 
post-processing step of deconvolving a final master image hypothesis to 
substantially correct a modulation transfer function describing a remainder of a 
distortion generated by the imaging sensor. 

5. The method of claim t s wherein said correction to the current 
master image hypothesis includes combining the correction images by deriving 
a weighted average of values of corresponding pixels in said correction images, 
the weight of each pixel in each correction image being calculated as a function 
of a distance as measured in the sampled image between: (i) a point in the 
sampled image to which the pixel in the correction image is mapped by the 
transformation operator; and (ii) at least one pixel centroid proximal to said 
point, 

6. The method of claim 5, wherein said function of a distance is 
derived from distortion operator BL 

7. A method for iterative derivation of a master image from a 
plurality of sampled images of non- identical at' least ' partially' overlapping, 
regions of a scene, the master image having an output resolution greater than a 
maximum resolution of the sampled images, the method comprising: 

(a) for each sampled image, defining: 
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(i) a transformation operator F mapping positions within the 
master image to corresponding positions in the sampled 
image, 

(ii) a distortion operator H simulating only a first portion of a 
distortion generated by an imaging sensor from which the 
sampled image was generated, and 

(iii) a sampling operator D for reducing an image from the 
output resolution to the resolution of the sampled image; 

(b) for each sampled image, applying said transformation operator, 
said distortion operator and said sampling operator to a current 
master image hypothesis so as to generate a predicted image, and 
calculating a difference image having pixel values corresponding 
to the difference in corresponding pixel values between the 
sampled image and the predicted image; 

(c) performing back-projection of each of said difference images to 
generate a correction image for the current master image 
hypothesis; 

(d) employing said correction images to perform a correction to the 
current master image hypothesis to generate a new master image 
hypothesis; and 

(e) after performing steps (b) through (d) at least once, 
deconvoluting a final master image hypothesis to substantially 
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correct a remaining portion of a distortion generated by the 
imaging sensor from which the sampled image was generated. 

8. The method of claim 7, wherein said correction to the current 
master image hypothesis includes combining the correction images by deriving 
a weighted average of values of corresponding pixels in said correction images, 
the weight of each pixel in each correction image being calculated as a function 
of a distance as measured in the sampled image between: (i) a point in the 
sampled image to which the pixel in the correction image is mapped by the 
transformation operator; and (ii) at least one pixel centroid proximal to said 
point. 

9- The method of claim 8, wherein said function of a distance is 
derived from distortion operator BL 

10, The method of claim 7, wherein said back projection includes 
employing an operator H bp corresponding to a pseudo-inverse of distortion 
operator H, wherein H bp approximates to an inverse of H at spatial frequencies 
below a given value and approaches zero at spatial frequencies above said 
given value. 
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11. The method of claim 7, wherein said back projection includes 
employing an operator H bp corresponding to a pseudo-inverse of distortion 
operator H, wherein H bp is chosen to substantially satisfy the condition: 



wherein: 



I is the unit operator for an image of the output resolution; 

D is a sampling operator for reducing an image from the output 

resolution to the resolution of an input image; and 
D l is an inflation operator for expanding an image from the resolution of 

the input image to the output resolution. 
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UPDATE OF REFERENCES 

-PCT/IL05/0001 16: International Search Report is enclosed. 
Additional applications were filed in: US, EP, CA, AU - no citations yet. 

CD including full copies of references cited is enclosed. 
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principle or theory un deriving the invention 
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Israeli Patent Office 



25 August 2009 

Notification of deficiencies in Israeli application 160152 

In correspondence with rule 41 of the Patent Rules, we hereby inform you that deficiencies 
have been found in the above-mentioned application. 

A response must be submitted within 4 months, extensions of time requested before the 4 
months* deadline at the cost of 59 NIS per month. 

1. The independent claims 1 and 7 ( version 2 ) are drawn to separate invention, 
contrary to clause 8 of the Patent law. 

Invention 1 - drawn to super division with backwards projection which approximate at low 
frequencies the opposite action of distortion which results from the imaging sensor and 
nears zero in remaining frequencies. 

Invention 2- drawn to super division without averaging at least part of the distortion 
resulting from the imaging sensor into a repetitive calculation with a compensation that is 
received after the processing step. 

Under Sec 24 (2) of the Law. You are required to divide the application. Please note that 
under rule 51 and Registrar Notice 62, the divisional applications should be filed within 4 
months or request for extensions should be submitted.. 

Sharon Perry, Examiner 



IN THE ISRAEL PATENT OFFICE 



27 December 2009 
Our File: 26/604 

Memorandum re: Super-Resolution Image 
Processing 

Israel Application No.160152 

Filed: 1 February 2004 

By: Rafael Advanced Defense Systems 

Ltd. 



Dear Sirs, 

This is a response to the Examiner's letter of 25 August 2009, which is being 
filed on or before 25 December and for which no extension fees are due. 

As per the examiner's request, claims 7-11 are being cancelled from this 
application, and are being submitted in parallel as part of a divisional application. 

Please find enclosed replacement pages 3-7 of the specification, in duplicate, 
which amend the summary to conform to the claims, as required by Hozer Rasham 
P23. 

Please also find enclosed replacement page 21, in duplicate, to which a 
sentence was added to ensure that the specification corresponds with the claims. 

Finally, please also please find enclosed replacement page 24 which is now 
the last page of the claims and to which a signature has been added. Pages 25-26 have 
been deleted. 



Respectfully submitted, 

Dr. Mark M. Friedman 
Advocate, Patent Attorney 
Moshe Aviv Towers, 54 th Floor 

7 Jabotinski Street 

Ramat Gan 52520 



160152/2 

partially overlapping, regions of a scene. The present invention also 
provides a technique, useful both alone and as part of the iterative derivation, 
for combining a plurality of input images into a single higher resolution output 
image. 

5 

According to the teachings of the present invention there is provided, a 
method for iterative derivation of a master image from a plurality of sampled 
images of non-identical, at least partially overlapping, regions of a scene, the 
master image having an output resolution greater than a maximum resolution of 

10 the sampled images, the method comprising: (a) for each sampled image, 
defining: (i) a transformation operator F mapping positions within the master 
image to corresponding positions in the sampled image, (ii) a distortion 
operator H simulating a modulation transfer function associated with an 
imaging sensor from which the sampled image was generated, and (iii) a 

15 sampling operator D for reducing an image from the output resolution to the 
resolution of the sampled image; (b) for each sampled image, applying the 
transformation operator, the distortion operator and the sampling operator to a 
current master image hypothesis so as to generate a predicted image, and 
calculating a difference image having pixel values corresponding to the 

20 difference in corresponding pixel values between the sampled image and the 
predicted image; (c) performing back-projection of each of the difference 
images to generate a correction image for the current master image hypothesis; 
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and (d) employing the correction images to perform a correction to the current 
master image hypothesis to generate a new master image hypothesis, 
wherein said back projection includes employing an operator H bp 
corresponding to a pseudo-inverse of distortion operator H, wherein H p 
approximates to an inverse of H at spatial frequencies below a given value and 
approaches zero at spatial frequencies above said given value. 



According to a further first aspect of the present invention, The method 
of claim Error! Reference source not found., wherein H bp is chosen to 
substantially satisfy the condition: 

H bp x D l x D x H = I 

wherein: 

I is the unit operator for an image of the output resolution; 

D is a sampling operator for reducing an image from the output 

resolution to the resolution of an input image; and 
D 1 is an inflation operator for expanding an image from the resolution of 
the input image . 



4 
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According to a further aspect of the present invention wherein distortion 
operator H corresponds to a combination of a modulation transfer function 
resulting from an optical system of the imaging sensor and a modulation transfer 
function resulting from a distortion generated by a sensor element array of the 
5 imaging sensor. 



According to a further aspect of the present invention, wherein distortion 

t 

operator H corresponds to a modulation transfer function describing only a first 
portion of a distortion generated by the imaging sensor, the method further 
lOcomprising a post-processing step of deconvoluting a final master image 
hypothesis to substantially correct a modulation transfer function describing a 
remainder of a distortion generated by the imaging sensor. 



According to a further aspect of the present invention wherein said 
15 correction to the current master image hypothesis includes combining the 
correction images by deriving a weighted average of values of corresponding 
pixels in said correction images, the weight of each pixel in each correction image 
being calculated as a function of a distance as measured in the sampled image 
between: (i) a point in the sampled image to which the pixel in the correction 



5 
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Image is mapped by the transformation operator; and (ii) at least one pixel centroid 
proximal to said point. 

According to a further aspect of the present invention, wherein said 
function of a distance is derived from distortion operator H. 

5 



10 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with 
reference to the accompanying drawings, wherein: 
5 FIG. 1 is a flow diagram of a method according to the teachings of the 

present invention for iterative derivation of a high-resolution master Image 
from lower resolution sampled images; 

FIG. 2 is a schematic illustration of the iterative process performed by 
the method of Figure 1; 
10 FIGS. 3 A and 3B are schematic representations of the sensor array of an 

imaging sensor and of an associated modulation (distortion) of imaging 
sensitivity, respectively; 



7 
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from the optical system MTK This decompilation may be performed by 
any of a number of known deconvolution techniques. The deconvolution process 
of the post processing is thus rendered very much more efficient since it Is 
performed only once on the final image rather than on each sampled image. 

It will be appreciated that the above descriptions are intended only to serve 
as examples, and that many other embodiments are possible within the scope of 
the present invention as defined in the appended claims. 



Material that exceeds the scope of the claims does not constitute part of 
the claimed invention* 
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4. The method of claim 1, wherein distortion operator H corresponds to a 
modulation transfer function describing only a first portion of a distortion 
generated by the imaging sensor, the method further comprising a post- 
processing step of deconvo luting a final master image hypothesis to 
substantially correct a modulation transfer function describing a remainder of 
a distortion generated by the imaging sensor . 

5 The method of claim 1 ? wherein said correction to the current master 
image hypothesis includes combining the correction images by deriving a 
weighted average of values of corresponding pixels in said correction images, 
the weight of each pixel in each correction, image being calculated as a 
function of a distance as measured in the sampled image between: (i) a point 
in the sampled image to which the pixel in the correction image is mapped by 
the transformation operator; and (ii) at least one pixel centroid proximal to 
said point. 

6. The method of claim 5, wherein said function of a distance is derived 
from distortion operator HL 



Dr. Mark Friedman Ltd. 
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4, The method of claim I , wherein distortion operator H corresponds to a 
modulation transfer function describing only a first portion of a distortion 
generated by the imaging sensor, the method further comprising a post- 
processing step of deconvolving a final master image hypothesis to 
substantially correct a modulation transfer Sanction describing a remainder of 
a distortion generated by the imaging sensor. 

5 The method of claim 1, wherein said correction to the current master 
image hypothesis includes combining the correction images by deriving a 
weighted average of values of corresponding pixels in said correction images, 
the weight of each pixel in each correction image being calculated as a 
function of a distance as measured in the sampled image between: (i) a point 
in the sampled image to which the pixel in the correction image is mapped by 
the transformation operator; and (ii) at least one pixel centroid proximal to 
said point. 

6, The method of claim 5, wherein said function of a distance is derived 
from distortion operator H. 
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WHAT IS CLAIMED IS: 



1. 



A method for iterative derivation of a master image from a 



plurality of sampled images of non-identical, at least partially overlapping, 



regions of a scene, the master image having an output resolution greater than a 
maximum resolution of the sampled images, the method comprising: 

(a) for each sampled image, defining: 

(i) a transformation operator F mapping positions within the 
master image to corresponding positions in the sampled 
image, 

(ii) a distortion operator H simulating a distortion associated 
with an imaging sensor from which the sampled image 
was generated, and 

(Hi) a sampling operator D for reducing an image from the 
output resolution to the resolution of the sampled image; 

(b) for each sampled image, applying said transformation operator, 
said distortion operator and said sampling operator to a current 
master image hypothesis so as to generate a predicted image, and 
calculating a difference image having pixel values corresponding 
to the difference in corresponding pixel values between the 
sampled image and the predicted image; 

■ 

(c) performing back-projection of each of said difference images to 
generate a correction image for the current master image 
hypothesis; and 
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(d) employing said correction images to perform a correction to the 
current master image hypothesis to generate a new master image 
hypothesis, 

i 

wherein said back projection includes employing an operator H bp 
corresponding to a pseudo-inverse of distortion operator H, wherein H bp 
approximates to an inverse of H at spatial frequencies below a given value and 
approaches zero at spatial frequencies above said given value. 

■ 

2. The method of claim 1, wherein H bp is chosen to substantially 
satisfy the condition: 

H bp x D £ x D x H = I 

wherein: 

I is the unit operator for an image of the output resolution; 

D is a sampling operator for reducing an image from the output 

resolution to the resolution of an input image; and 
D l is an inflation operator for expanding an image from the resolution of 

the input image to the output resolution, 

3. The method of claim 1, wherein distortion operator H 
corresponds to a combination of a modulation transfer function resulting from 
an optical system of the imaging sensor and a modulation transfer function 
resulting from a distortion generated by a sensor element array of the imaging 
sensor. 
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4. The method of claim 1 5 wherein distortion operator H corresponds to a 
modulation transfer function describing only a first portion of a distortion 
generated by the imaging sensor, the method further comprising a post- 
processing step ■ of deconvoiuting a final master image hypothesis to 
substantially correct a modulation transfer function describing a remainder of 
a distortion generated by the imaging sensor. 

5 The method of claim 1, wherein said correction to the current master 

■ 

image hypothesis includes combining the correction images by deriving a 
weighted average of values of corresponding pixels in said correction images, 
the weight of each pixel in each correction image being calculated as a 
function of a distance as measured in the sampled image between: (i) a point 
in the sampled image to which the pixel in the correction image is mapped by 
the transformation operator; and (ii) at least one pixel centroid proximal to 
said point. 

6, The method of claim 5, wherein said function of a distance is derived 
from distortion operator BL 
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